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Abstract 

We consider the effects of a velocity-independent friction force on 
cantilever damping. It is shown that this dissipation mechanism causes 
nonlinear effects in the cantilever vibrations. The size of the nonlinear- 
ity increases with decreasing cantilever velocity. Our analysis makes 
it possible to understand Stipe's et al. pQ experiments where an am- 
plitude dependence of the cantilever eigenfrequency and anomalous 
dissipation was observed only at small amplitudes. 



1 Introduction 

The quality and potential practical opportunities for micromechanical de- 
vices strongly depend on their dissipative characteristics. Numerous authors 
have studied the damping of cantilever vibrations both theoretically and ex- 
perimentally. There are at least two obvious motivations for these studies. 
First, the amount of dissipation determines the sensitivity of cantilever-based 



"Corresponding author: gpb@lanl.gov 



1 



devices. Hence, reducing dissipation improves the technical characteristics of 
the devices. The second is based on the possibility of probing the dissipative 
forces which arise due to cantilever-sample interaction in order to obtain a 
useful information about the sample. The last idea was a basic one for orig- 
inating a new direction for atomic force microscopy: noncontact dissipation 
force microscopy. (See references |2]-[7j). 

Many attempts to understand the dissipation mechanisms in cantilevers 
were undertaken. The corresponding theoretical models, which could de- 
scribe cantilever damping, are presented in references [8] and [9]. At the 
same time, the fundamental nature of damping is still not clear, especially 
in the case of noncontact cantilever vibrations. (See, for example, discus- 
sions in references [I], [2], [7], an d [13] )• Fluctuations of the van der Waals 
force ( "vacuum friction" ) or Joule losses are not responsible for the dissipa- 
tion measured experimentally. The corresponding friction effects have been 
calculated in [10]- [13] to be many orders of magnitude smaller than those 
measured in p]. There is not only quantitative disagreement, but also differ- 
ent qualitative behavior. The point is that both theories predict very strong 
dependence of damping on the cantilever-sample separation, and this is not 
consistent with the experimental observations. 

A more complex and more specific model, involving adsorbed particles, 
was suggested in [TJ]. The authors connect the large long-range noncontact 
friction with the electromagnetic interaction of moving charges, induced on 
the surface of the tip by the bias voltage, with acoustic vibrations in an 
adsorbate layer on the surface of the sample. 

An internal quantum friction due to the presence of a finite number of 
two- level systems is analyzed in Refs. [15] and [16] . Their analysis is related 
to millikelvin temperatures which is beyond experimental conditions of PQ. 

The experimental dependence of the friction force, F, on the bias voltage, 
V, applied between the tip and a sample, has been found. It is the quadratic 
dependence, F ~ V 2 (see, for example, references [1] and [7J) that is really 
typical of the Joule mechanism for energy losses. At the same time, it should 
be emphasized that not only Joule losses, but also a tip-sample attractive 
forces, caused by electrostatic interactions between closely spaced charged 
surfaces, behave as V 2 . Taking into account the fact that the Joule mecha- 
nism is not effective (at least for good conductors), the increased dissipation 
may be attributed to some effects of the attractive forces. In this connection, 
it should be noted that the increase of dissipation with decreasing tip-sample 
spacing was also obtained in the absence of the electrostatic attraction (for 
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zero- value bias- volt age). The calculations of [T3] show that at small separa- 
tions the van der Waals (Casimir) attractive force may become so strong as 
to be able to modify the cantilever eigenfrequencies. 

A very specific attraction of the metal tip to dielectric samples occurs 
as well. Although overall the dielectric sample is electrically neutral, it may 
contain localized electric charges, which interact with the corresponding "im- 
age" charges in the tip. Samples of fused silica were used by authors of [T] to 
study the effect of charged centers on the cantilever vibrations. (Reference 
[T3] contains the corresponding theoretical calculations.) The number of cen- 
ters has been varied by means of irradiation with 7 rays. An enhancement 
of dissipation occurs when the concentration of charged defects increases. 

The increasing attractive forces in the above systems are accompanied 
by the growth of dissipation. Although this force may be of a different 
physical nature, it always results in an increase in cantilever dissipation. 
Therefore any consistent theory of dissipation should take into account this 
circumstance. 

Concluding the introduction, we highlight a statistical analysis of various 
micromechanical systems [5] that shows that the quality factor (inverse of 
dissipation) usually scales linearly with size. There is a general tendency of 
dissipation to increase with the surface/volume ratio, indicating the relevance 
of the cantilever surface layers to dissipation. 



2 The friction force 

In general, the cantilever motion is described by the oscillator model. The 
corresponding equation of motion is given by 

where X = X(t) is the displacement of the cantilever tip; m and k are the 
cantilever effective mass and spring constants, respectively. The term is 
the friction force, which is usually assumed to be proportional to the velocity 
of the cantilever motion. Hence, this force tends to zero when ^ — > 0. 
The experiments in [TJ were performed for the cantilever frequency ujq = 
(k/m) 1 / 2 « 24 x 10 3 s _1 and vibration amplitudes of the order of lOnm or less. 
This corresponds to the characteristic value of the velocity equal to 0.02cm/ s. 
This slow cantilever motion, in practice, excludes dissipation mechanisms due 
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to both Joule losses and fluctuations of the van der Waals forces. Therefore, it 
seems to be quite reasonable to consider a velocity-independent friction force, 
F, as the one responsible for dissipation in this case. The conventional term 
"velocity-independent force" means a force which does not depend on the 
absolute value of the velocity but always acts in the direction opposite to the 
tip velocity. Thus, F reverses its sign after each half-period of vibrations. 
This force occurs in mechanical systems in which one solid surface is in 
contact with another and slides along it. 

The friction force is proportional to the effective contact area, which, 
depends linearly on the compressing force (the Coulomb friction law) . A very 
exotic mechanism of velocity-independent friction due to photon exchange 
between surfaces which are in close proximity to one another (but not in a 
contact) was predicted in Ref. [TT] . 

For the case of a cantilever pQ whose body has a complex multilayer 
structure, an internal friction force between the different layers may arise. 
The body is a platinum-coated single- crystal silicon beam. There is also 
a lnm titanium layer between the lOnm platinum "coat" and the silicon. 
Besides that, the tip was coated with 200nm of gold by evaporation. It 
was established experimentally that the dissipation of the platinum-coated 
single-crystal silicon cantilever is higher by a factor of 6 than that of an 
uncoated cantilever. Therefore, a muli-layer structure is essential in forming 
the overall dissipation. It seems to be reasonable to consider a tip-sample 
stretching force to be responsible for increasing the contact interlayer area 
and, consequently, the frictional force, both of which linearly depend on the 
tip-sample attraction. 

The cantilever motion is described by the equation 

+ = F, (2) 

where F > for ^ < and F < for ^ > 0. Formally, Eq. [2] may be 
rewritten in the form given by Eq. [TJ with the coefficient V depending on the 
velocity: 

The oscillator equation with V given by Eq. [3] is no longer a linear equation. 
Nevertheless, it can be easily solved for each time interval between any two 
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Figure 1: Ringdown oscillations in the case of a velocity-independent force. 
Points, where X = 0, are slightly displaced from the values n + ^ which 
would correspond to harmonic oscillations. 



successive turning points. The general solution is given by 

X(t) = (-If + [-X + ^(2n + l)]cos(u t), (4) 

where n is the integer part of Uot/n and — X < is the initial position of 
the tip. Eq. H] describes the decaying motion of the cantilever. In the range 
of small dissipation (the energy loss per period is small compared with the 
total energy e) and small values of time, Eq. @] may be approximated by a 
simplified formula 

X{t) w -X (l - at)cos(uj t) w -X e~ at cos(u t), (5) 

where 

2\F\ 

a = a(X ) = — ■ (6) 

irmuj X 

The coefficient a describes the nonlinear damping of the oscillations, and it 
depends strongly on the amplitude Xq. As we see, the nonlinearity is the 
most pronounced for small amplitudes as was observed experimentally [1] (at 
amplitudes of the order of or less than lOnm). It was noticed in Ref. [T] that 
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this effect vanishes if either the amplitude of oscillation or the tip-sample 
separation increases. 

With the known oscillator trajectory described by Eq. it is easy to 
estimate the value of nonlinear shift of the oscillation frequency. In what 
follows we consider the so called "ringdown" regime of oscillations, which 
takes place when a drive circuit is abruptly grounded and the cantilever rings 
down until thermal equilibrium is established. The value of the cantilever 
frequency may be obtained by measuring the time interval t 2 — t\ between 
two successive crossing points where the displacement X is equal to zero. For 
our model, these points do not match with the values it/2ujq and 3n/2u as 
would occur for harmonic oscillations (see Fig. 1), but are slightly displaced 
to greater values of t. Considering these displacements as small quantities, 
we can easily obtain the difference t 2 — t\. 

Let us introduce the notation, co> ti = § + 5\ and Uot 2 = + S 2 . Then 
we obtain from Eq. (TJJ two equations for 5\ and 5 2 . They are 

\F\ ( \F\\ 
— — h — X H — — sindi = Oand 



k V k 

^ + ( - X + 3^) sin5 2 = 0. (7) 

It follows from Eqs. © that ^ « g(l + g) and S 2 « g(l + 3g). 
Then we have 

ujq Uok 2 X l 

The corresponding value of the frequency is given by 



t 2 -t! = — + 12v2 . (8) 
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2F 2 



where the energy of vibrations, e, is equal to kX^/2. 

The experimentally observed decrease of the frequency is illustrated in 
Fig. 2d of reference [1] . The decrease in the oscillation frequency takes place 
at the beginning of the transition from high- to small-amplitude regimes of 
oscillations. It would be interesting to study experimentally the dependence 
of the frequency shift on the tip-sample attraction, F attr . For the case of 
the Coulomb law for the friction force, i.e. when |F| = c + C\F attT (c i are 
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costants), the frequency shift is given by the expression 




) 
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Co + CiF a ttr 



(10) 



When the attraction force is induced by the bias voltage, V, F attr ~ V 2 . 
When it is due to localized charges in the dielectric sample, F attr ~ n ch , 
where n c h is the concentration of charges. In both cases the model proposed 
here can be examined experimentally. 

There is an alternative method of calculation of the frequency shift. It 
is based on the asymptotic theory of nonlinear oscillations developed by 
Bogolubov and Mitropolskii in [18J. The direct application of this method to 
our specific case of a velocity-independent friction force results in a frequency 
shift which is slightly smaller (by factor of 3/n) than that given by Eq. 
The details are in Appendix A. 

So far our analysis assumed a constant friction force. This force may 
be considered as originating from many discrete events of energy losses. For 
brevity, we will use the term "kick" for each such event. The duration of each 
kick is assumed to be very short compared with the period of vibrations, and 
the number of kicks per period is large. After averaging over a time interval 
much shorter than the period of oscillations, but large enough to include 
many kicks, the effect of kicks reduces to a constant friction force. It is 
evident that the description based on a continuous friction force is applicable 
to high-amplitude oscillations only. The discrete nature of the dissipation 
may reveal itself at small amplitudes when the mean "free path" I is longer 
than or of the same order as the amplitude. The correspondence of the 
two models to each other imposes a relation between the parameters of both 
\F\ = Ae/Z, where Ae is the energy loss in the course of each kick (Ae << e). 

In what follows, we consider the case of a low-probability kick per oscil- 
lation period. This is possible when X << I. If the kick occurs during the 
motion away from (towards) the center, the time of returning to the central 
position (X = 0) is decreased (increased). Hence, the period of oscillations 
is a fluctuating quantity. It can be shown that these fluctuations result in 
a positive frequency shift that is in contrast to the case of large amplitudes, 
X » I. The following simple analysis makes it possible to obtain the value 
of the shift in an explicit form. As before, we proceed from the equation of 
motion for a harmonic oscillator 



X(t) = X cos(uot) 



(11) 
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which describes the vibrations in the absence of dissipation. Initially, the 
oscillator is in the position X(t) = X . Let us assume that the kick occurs at 
time t\. Then X(t\) = X\ — X cos(uJoti). The value of t\ can be expressed 
via X\ as 

t x = — cos — 

\X , 

After this event, the oscillator continues its motion with a modified amplitude 
C, which can be obtained from energy conservation. It is given by 



C = ^Xq 2 - 2Ae/k. (12) 

For simplicity, we will consider that Ae << e. 
The modified equation of motion is given by 

X(t) = Ccos{u t + <p), (13) 

where if is the phase variation introduced by the kick. It is convenient to 
introduce the new time variable, r = t + (p/uo. It can be easily seen that r 
varies from "Cos -1 ^ 1 (the instant of kick) to ^ (the instant of crossing the 
x axes). Then the overall time required to get to the central position X = 
is given by 

7T 1 7T , 

h + - cos — = — + r + , 

2uj uj C 2(jj 

where 

It follows from Eq. [TD that the time required to reach the center has 
increased due to the kick. In contrast, the transit time decreases when the 
oscillator moves away from the center. The corresponding time is given by 

+ where 

r~ = — ( sin" 1 - - — sin' 1 — - ). (15) 

As we see the variation of the transit time depends essentially on the position, 
Xi, of the oscillator at the instant of kick. We assume that all realizations 
of Xi have equal probabilities. This assumption will make it possible to take 
into account the contribution of the transit-time fluctuations to the frequency 
shift. 
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As before, our following analysis concerns the case of the ringdown ex- 
periment. The time required to undertake N full cycles can be expressd 
as 



t 



u 



A), 



(16) 



where 



A 



2nN 



£ t~ (Xi)rn + Yl T+ ( x j) n j 



The indices i and j denote the quarter-periods with away-from-center and 
towards-center motion, respectively. The random variables equal to 

1 or depending on whether or not a kick occurs during the (i,j)th quarter- 
period. The average value of n« is given by (rii) = ^P- « 1, implying 
a negligible probability of two kicks during any given quarter-period. The 
period of oscillations, T, is defined as T = t/N. Hence, the frequency of 
oscillations is given by 



2nN/t = u (l + A)- 1 « u (l - A + A 2 - ...) 



(17) 



After statistical averaging of lu, the linear term (A) vanishes. This follows 
from the definitions of r^ 1 and the trigonometry relation sin~ 1 {x) = | — 
cos~ 1 (x). The main contribution to the frequency shift is from the term 
(A 2 ). Thus, we have 



(u) - uj 



7C 2 N 



(18) 



In deriving Eq. EDI we use the following relations: 



i nt = Hi 



(mnj) 



{rii){nj) when i ^ j. Also, the integer N cannot be very large because the 
energy dissipation during the observation time is assumed to be small. 

Considering the losses as statistically independent events which have 
equal probabilities at any value of X within the interval [— C, +C] (the in- 
terval where a kick is possible), we can calculate the average value of [t~] 2 
as 



■ ^( x 



sin 



-i 



X 
X~ 



This integral can be calculated analytically considering 4r << 1- 



Appendix B). Finally, we get 

(uj) -uj r 



2c 



\F\Ae 



AV2tt 2 N 



me 3 / 2 \AeJ 



(19) 
(See 

(20) 
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Comparing the amplitude dependence of the frequency shift in two limit- 
ing cases given by Eqs. [H]and[2UJ we may conclude that the function u(X ) 
has a minimum at intermediate amplitudes as shown in Fig. 2. This curve 




Figure 2: Dependence of frequency shift on the amplitude of oscillations Xq. 
is similar to the experimental one displayed in Fig. 2d of reference pp. 

3 Conclusion 

We have described phenomenologically a possible mechanism that produces a 
small-amplitude nonlinearity for cantilever vibrations. It is assumed that the 
nonlinearity is due to velocity-independent friction force which, at small am- 
plitudes, reduces to a sequence of random kicks retarding cantilever motion. 
Our analysis predicts both the increase of the decrement of oscillations and 
nonmonotonic behavior of the frequency shift when the amplitude decreases. 
It is shown that the last effect is due to the crossover from the amplitude- 
independent friction force to discrete events of the cantilever energy losses. 
This is in a qualitative agreement with the cantilever vibrations observed 
experimentally. 
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Appendix A 

A regular method for solution of the equations 

X(t)+u 2 X = ef(X,X), (21) 

where e is a dimensionless quantity (e — ► 0), is developed in [18]. In our case, 
the term on the right-hand side is given by 

ef(X,X) = --£-. (22) 
m \X\ 

According to [18] , the solution of Eq. [21] can be represented as an expansion 
in powers of e, 

X{t) = acosip + eui(a, <p) + e 2 w 2 (a, <p) + 0(e 3 ), 

a = eAi(a) + e 2 A 2 (a) + 0(e 3 ), 

<p = LO + eB 1 (a) + e 2 B 2 (a) + O(e 3 ). (23) 

To obtain the decrerment and frequency shift of oscillations, it is sufficient 
to calculate only terms up to the second order in e 2 . For the case of the 
dissipation force given by Eq. [221 these quantities are determined by 

Ai(a) = A, 

B x (a) =A 2 (a) = 0, 

*w = -s#-s*?S^' (24) 



where 



4|/| . fnn 
-sin\ 



nn V 2 
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The quantity A\ determines the decay of oscillations. Setting the initial 
condition for the amplitude a corresponding to the solution given by Eq. HJ 
a(t = 0) = — Xq, we obtain from Eqs. [231 and [2H 



a ^--H l -^xV- (25) 

This coincides with the results given by Eqs. [5] and [61 

The value of the renormalized frequency u(a) = (fir a) can be obtained 
explicitly after summing up in the expression for -62(0). This sum is equal 
to F 2 /4. Then the frequency of oscillations is given by 

/ QF 2 \ 
" = MX$)- ( 26 ) 

The second term in the brackets is very close to that in Eq. [91 Their ratio 
is equal to 3/tt which is close to unity. Thus both approaches give almost 
identical results. 



Appendix B 



1 r c 



({r-} 2 ) = cJ Q dX 



l( X\ ,(X 



-\ 2 







C) \x , 

dx[sin~ l x — sm~ 1 x(l — 5)} 2 , (27) 
where 5 = Ae/(2e) « 1. Let us set 

sin~ 1 x(l — 5) = y + 5' , 
where \S'\ << 1. Then we have siny = x and 



8' ss Vx- 2 - 1 - Vx- 2 -1 + 25. 
After substituting x = (1 + z)' 1 ^ 2 , the last integral in (1271) is reduced to 

r°° dz 1 



(l + zf/ 2 z + 5 + Vz 2 + 2z8' 



(28) 



12 



One cannot set S = in the integrand of Eq. (1281) because of logarithmic 
divergence of the integral. 

Let us denote the last integral as K. It is convenient to divide the range of 
integration in two parts and express the whole integral as a sum K = K1+K2, 
where 



The integrand in K\ can be approximated by [z + 5 + \J z 1 + 2zS] 1 , and it 
can be easily seen that K\ ~ (l/4)|/n<5|. 

r 1 _1 

The integrand in K2 can be approximated by (1 + z) 3 ^ 2 2z . For small 

values of S, K 2 « [1/A)\ln5\. Then, K w {l/2)\ln5\. Finally, using Eqs. 
and (128p . we have 

<M 2 > - ^<«|. (») 
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